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Abstract. We present an analysis of the location of inter- 
acting E/SO galaxies on the Fundamental Plane (FP). Us- 
ing the NEMO package, we performed N-body simulations 
of close encounters and mergers of two self-gravitating 
spherical galaxies. Two models for encounting galaxies 
- Plummer's and Hernquist's - were used. The number 
of particles ranged from 20000 to 50000 per galaxy. The 
changes of central density, half-mass radius and central 
velocity dispersion were analysed. It was found that close 
encounters between galaxies alter noticeably the above pa- 
rameters within a very short time interval (10 7 -10 8 years) 
just before the final merger. The amplitudes of parame- 
ter changes strongly depend on the initial mass concen- 
tration of a model. In some experiments we considered 
j an encounter of two spherical galaxies with dark matter 

■ components. The effect of parameter changes was less pro- 
nounced than for the experiments without a dark halo. 

• The results of the simulations were used to discuss the FP 
I ' for interacting early-type galaxies. 
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1. Introduction 

The Fundamental Plane (FP) defines one of the most im- 
portant relationships for early-type galaxies (Djorgovski & 
Davis 1987, Dressier et al. 1987). The FP combines surface 
photometry parameters (R e - effective radius and \i e - 
effective surface brightness or (/i) e - mean surface bright- 
ness within R e ) with spectroscopy characteristics (line-of- 
sight central velocity dispersion cr ). The measured values 
of R ei /i e and cr for a sample of E and SO galaxies do 
not fill this 3-paramctcr space entirely but rather a thin 
plane within it (with scatter ~ 0.1 dex). The FP can be 
projected onto any two axes. Examples of these projec- 
tions are the Kormendy relationship (n e -\ogR e ) , and the 
Faber-Jackson relationship between luminosity and veloc- 
ity dispersion. 

The interpretation of the FP is still a matter of dis- 
cussion. The most popular point of view is that the FP is 
simply a consequence of the Virial theorem and the fact 
that E/S0 galaxies have similar mass-to-luminosity ratios 
(M/L) and homologous structure at a given luminosity. 
If elliptical galaxies had perfectly homologous structural 
and kincmatical properties, then the Virial theorem would 
predict the FP to have the form R e oc er^E^M/L) -1 , 
where E e = 10~°- 4 ^ e . Under the assumption of a con- 
stant M/L, the observed FP should then be logR oc 
21ogcr + 0.4(^) c . In practice, however, the logcr and [i e 
coefficients have been determined to lie in the range 1.2- 
1.6, and 0.30-0.35, respectively. The difference between 
expected and observed parameters of the FP is generally 
refered to as the "tilt" of the FP. This tilt is usually viewed 
as the product of systematic variations of Mj L along the 
FP (M/L oc L@). Among the possible explanations that 
have been considered for the tilt are: variations in the stel- 
lar population along the FP, changes in the dark matter 
content, and deviations from homology (Hjorth & Madsen 
1995, Capelato et al. 1995, Ciotti et al. 1996). 

The very existence of the FP has profound implications 
for the processes that led to the formation and evolution 
of elliptical galaxies. In particular, its tightness provides 
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a very strong constraint on the evolutionary history of 
those systems. The FP can be used to study the evolu- 
tion of galaxies as a function of redshift, to derive rcdshift 
independent distance estimates, etc. 

However, we must be cautious when using the FP to 
obtain distances to elliptical galaxies without regard for 
their environment. The significant environmental effects 
on the elliptical galaxy correlations prevent their utility as 
accurate distance indicators. Hence the question of how a 
galaxy should move in the space of (R e , fi e , <r ) as a result 
of interaction with other galaxies is of a great importance. 

Zepf & Whitmore (1992) investigated the environmen- 
tal dependence of the FP of elliptical galaxies by studying 
a sample of elliptical galaxies in compact groups. The cen- 
tral velocity dispersions of ellipticals in compact groups 
appear to be about 20 % lower than those in other environ- 
ments (galaxies in clusters and in the general field). De la 
Rosa et al.(2001) re-examined this topic with a better and 
more homogeneous dataset. They found little or no sig- 
nificant difference between elliptical galaxies in compact 
groups and comparable galaxies in other environments. 

Levine & Aguilar (1996) carried out a few Monte Carlo 
simulations of the evolution of a number of elliptical galax- 
ies in the core of a very rich cluster. Under two very simple 
physical assumptions (virial equilibrium and retention of 
the functional form of the surface brightness and mass 
density profiles), any change in the galaxy caused by any 
interaction that satisfies these assumptions will be nar- 
rowly constrained within a plane very close to the observed 
FP of ellipticals, even for large changes in the global M/L 
ratio. The angle between this plane and the observed FP 
of ellipticals is less than 15°. So, if a galaxy is in or near 
the FP, it will remain so after an interaction. 

De Carvalho & Djorgovski (1992) found that elliptical 
galaxies in clusters were best fitted by a slightly different 
FP than galaxies in the field. However according to Pahre 
et al. (1998) this difference is naturally explained by er- 
rors in estimating dust extinction or estimating distances 
(and not by any intrinsic differences in the stellar popu- 
lations among elliptical galaxies that correlates with their 
environment). There is no clear evidence that either the 
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Fig. 1. The Faber- Jackson relation for 94 interacting E/S0 
galaxies. 



slope, zero-point or scatter of the FP are dependent on the 
spatial environment (see discussion in Pahre et al. 1998). 

The main purpose of the present work is to explore the 
possible influence of strong mutual low-velocity encounters 
of early-type galaxies on their general characteristics and 
on their location within the FP. 



2. Fundamental Plane for Interacting and 
Merging Galaxies 

In order to study the FP for intensively interacting early- 
type galaxies we have considered several pieces of obser- 
vational data. 

Fig.l presents the Faber- Jackson relation for E/S0 
galaxies belonging to 1) VV (Vorontsov-Velyaminov 1959) 
or Arp (Arp 1966) systems (28 galaxies, stars), 2) CPG 
(Karachcntsev 1972) or TURNER binary systems (Turner 
1976), and 3) Hickson compact groups (Hickson 1982) (cir- 
cles). The total number of galaxies is 94. 

We used <7o from 1) Hypercat, CRAL Observatoire 
de Lyon (www-obs.univlyonl.fr/hypercat), 2) Davoust & 
Considere (1995), 3) Simien & Prugniel (1997a), 4) Simicn 
& Prugniel (1997b) and apparent magnitudes raj from 
NED0. The raj have been corrected for galactic extinc- 
tion using NED (reddenings derived from HI and galaxy 
counts: Burstein & Heiles 1982). Absolute magnitudes Mb 
are computed assuming Hq = 75 km s" 1 Mpc -1 . When es- 
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timating galaxy distances, we assigned to each galaxy the 
heliocentric radial velocity adopted from NED and cor- 
rected for the Sun's motion using the standard relation. 

Dashed line in Fig.l shows the mean relation for 594 
E/SO galaxies according to McElroy (1995): Lb k CTq' 3 - 
The slope and the scatter (cr(logcro) = 0.15) of the data 
for interacting galaxies are the same as for normal E/SO. 

Triangles in Fig.l present the characteristics of 14 rem- 
nants of disk-disk mergers (Keel & Wu 1995). As one can 
see, the mergers follow the same relation as early- type 
galaxies but with small shift (~ 0. m 6) to brighter ab- 
solute luminosities. The dashed rhomb shows the mean 
value (±lcr) for 9 starbursting infrared-luminous galax- 
ies (Shier & Fischer 1998). All these galaxies are at some 
stage of merging. Starbursting mergers demonstrate sig- 
nificant brightening (~ 2 m ) in comparison with normal 
galaxies. One can expect that relics of disk-disk mergers 
and infrared-luminous mergers will evolve into normal el- 
lipticals at the FP. 

Fig. 2 compares the FP location of interacting el- 
lipticals/lenticulars belonging to VV and Arp systems 
(crosses) with relatively isolated galaxies (points) (the 
data in the rQ passband are from Djorgovski & Davies 
1987, in the B filter - Bender et al. 1992). Here we used 
H = 100 km s- 1 Mpc" 1 . The FP is defined as logR oc 
logfTo + 0.26(/i) c . The FP scatter for interacting galaxies 
(a(logRe) = 0.09 in B and cr(logR c ) = 0.14 in r G ) is the 
same as (or even smaller than) that of non-interacting ob- 
jects (a(logRo) = 0.19 in B and <r(logR c ) = 0.15 in rc)- 

One can conclude that the present observational data 
about the global characteristics of strongly interacting 
E/SO galaxies show that the close environment has no or 
almost negligible effect on the FP. The only difference we 
see for forming (or young) ellipticals is that they are signif- 
icantly brighter (for fixed <7o) due to superimposed bursts 
of star formation. 



3. Numerical Modeling of Interacting Ellipticals 

To check our empirical findings we have performed numer- 
ical modeling of close encounters of two ellipticals. 

We simulated the evolution of the self-gravitating 
spherical galaxies by using the NEMO package 

1995) 



( bttp: / /bima. astro. umd / nemo / 



Teuben 

is a freeware package designed to numerically 



This 
solve 



gravitational N-body problems. It consists of subroutines 
for specifying initial configurations of stellar-dynamical 
systems (including many standard models) and subrou- 
tines for simulating the evolution of these systems based 
on various numerical schemes. 

In our computations, we used N-body dynamics with 
a hierarchical tree algorithm and multipole expansion to 
compute the forces, as proposed by Barnes & Hut (1986). 
Multipole expansion is limited by the monopole term. A 
tolerance parameter ("opening angle") which controlled 
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Fig. 2. The FP for interacting E/SO galaxies. 



the force computation from distant particles was set as 
0.7. 

We considered only the behavior of the stellar com- 
ponent of galaxies (elliptical galaxies usually contain a 
small amount of gas). The number of particles used in 
our numerical simulations ranged from N tot — 20 000 to 
Ntot = 50 000 per galaxy. In this case we managed to sub- 
stantially suppress the effects of pair relaxation and to 
trace the evolution of merging galaxies on time scales up 
to t « 0.5 x 10 9 years. The timestep for the simulations 
used was as At = 10 6 years. 

The softening length e was taken as < 0.1 of the inter- 
particle distance of the homogoneous state. 

We considered two models for encounting galaxies. One 
of them is Plummcr's spherically-symmetric model 



$(r) 



p{r) = 



GM 



(r2+a 2 ;) l/ 2 

3M 



47T 
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where M is the galaxy mass and a p i is a scale length. 

In some experiments galaxies are modeled by the 
potential-density pair proposed by Hernquist (1990) for 
spherical galaxies 



*(r) = - 



GM 
r + a hq 



M 



(r) = — 

P{ ' 2^al q r{r + a hq f 



hq 

It is well known that both models are described by 
a distribution function in an analytical form and, in the 



absence of numerical errors and dynamical instabilities, 
remain time-stationary. 

Results are presented in the following system of units: 
gravitational constant G = 1, the galaxy mass M = 1, the 
total energy of a sphere E = -1/4 (E = -'SttGM 2 /64a p/ 
and a p i — 3n/16 for Plummer's sphere; E = —GM 2 /I2ah q 
and a) lq = 1/3 for Hernquist's sphere). Scaled to physical 
values appropriate for a typical elliptical galaxy, i.e. M = 
10 11 and half-mass radius r 1 / 2 — 3 kpc (r 1 / 2 ~ l-31a p / 
- for Plummer's sphere, r\/ 2 ~ 2Alah q - for Hernquist's 
sphere), units of distance, time and velocity are 3.73 kpc, 
10.3 Myr, 345.3 km s _1 respectively 

We specified the initial distance between two equal- 
mass galaxies as R — 37.3 kpc and chose the initial relative 
velocity in the range V = 77.3 — 103.6 km s _1 . 

As a result we had a close encounter with merging 
(the distance of the first closest approach was 5.2 kpc) 
and a distant encounter without merging (in this case the 
minimum galaxy separation was 10.3 kpc). 

Fig. 3 presents some "snapshots" of a close encounter, 
showing the initial condition (i = 0), the configuration 
near the first maximum overlap (t — 30), the configu- 
ration shortly before the final merger (i — 34), and a 
merger state (t = 40). A contour map of two interact- 
ing Plummer's spheres are plotted in Fig. 4. Some changes 
in morphology of model galaxies become noticeable only 
at the final stages of encounter (i = 30 — 34). There are 
real objects which are very similar to these interaction 
stages of model galaxies, for example NGC 1587/1588, 
NGC 7236/7237 (Borne & Hoessel 1988). 

For all models we followed the changes of central den- 
sity (N - number of particles within 0.2 distance units 
= 0.75 kpc), half-mass radius (r 1 / 2 - within 0.75 kpc) 
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Fig. 3. Close encounter of two identical Plummer's spheres. Times: 0,2,6,14,28,30,32,34,40 (left to right, top to bottom). 



and central velocity dispersion (uo)- During the distant 
encounter all these parameters were not altered. As to 
the close encounter, in this case the most drastic changes 
in parameters were observed at the moment shortly before 
the final merger (at A t « 5 « 50 x 10 6 yrs before merger) 
- Fig.5. 

The range of parameter changes appears to depend on 
the initial mass concentration of a model. More concen- 
trated Hernquist's spheres show less pronounced responses 
than does a less concentrated Plummer's model (Fig.5). 

Amplitudes of relative changes of the parameters for 
two models are: 

Plummer's sphere - 8r\/2 ~ 67%, 5<tq « 17%, 
Hernquist's sphere - 5ri/ 2 ~ 33%, <5oo ~ 6%. 
These amplitudes are comparable with the FP scatter 
(rms scatter ~0.1 - 0.15 dex). Therefore even the close 
low-velocity encounters cannot significantly enlarge dis- 
persion of the observed FP. Moreover the galaxies demon- 
strate large changes only within a very short time interval 
(~T0 7 -10 8 years) just before the final merger (short rela- 
tive to the evolution time of merging galaxies). 

It should be noted that Hernquist's sphere shows large 
parameter deviations from initial values (more than the 



values given above) at the very last stages of the encounter 
before merger (t = 45 — 49) - Fig.5. We do not discuss 
these moments because of the short duration of these pe- 
riods and the difficulty in distinguishing almost merged 
galaxies. 

We have also considered an encounter between two 
galaxies embedded in dark halos. Each galaxy was com- 
posed of a spherical component, represented by a Plum- 
mer's model, and a spherical Plummer's halo. The mass 
of dark matter was set equal to the visible mass within 
the sphere containing 95 % of the stellar mass. The com- 
bined potential of both components was used to solve the 
Jeans equations and to derive the initial velocity disper- 
sion (Hernquist 1993). Each sphere consisted of N tot — 
40 000 particles (N tot = 20 000 per stellar component and 
Ntot = 20 000 per halo). 

The addition of an extended dark halo surrounding 
each galaxy makes the effect of parameter changes less 
pronounced than for the experiments without a dark 
halo. Amplitudes of relative changes of the parameters 
are: 

Plummer's sphere - Sr 1 / 2 « 60%, Sao ~ 25%, 
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Fig. 4. Contour maps of two interacting Plummer's spheres. Times: 0,6,14,28,30,32,34,40 (left to right, top to bottom). 
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Fig. 5. Evolution of the characteristics of two model galaxies (solid and open rhombs) during a close encounter 
(-/Vj ot =50 000 per galaxy). Top three figures - Plummer's spheres, bottom - Hernquist's spheres. Crosses show param- 
eters after merging of two Plummer's spheres in one object. Evolution of the corresponding isolated spheres are shown 
by solid lines, t = corresponds to initial configuration, £=30, 32 - first pericenter passage, i=34,38 - the time of 
largest separation after the pericenter passage and t «40,49 - merging (first and second numbers refer to Plummer's 
and Hernquist's models correspondingly). 



Plummers's sphere with halo - 5r 1 / 2 « 30%, Sao ~ 15%. 



4. Conclusions 

Our numerical experiments and observational data have 
shown that the global parameters of early-type galaxies 
are rather stable to gravitational perturbation. The FP 
parameters do change during close encounters but within 
a very short time interval (~10 7 -10 8 years) just before 
the final merger. Furthermore, the amplitudes of these 
changes are comparable to the scatter of the observed FP 
of ellipticals. There is a very small probability that we 
observe the system at these interaction stages. Therefore, 
the FP is still a good method to derive distances, even 
to the clusters that include a large number of interacting 
galaxies (e.g. van Dokkum et al. 1999). 

Acknowledgements. We acknowledge support from the "Inte- 
gration" programme (A0145). We would like to thank the ref- 
eree, Dr. R. de Carvalho, for useful comments and suggestions. 



References 

Arp, H. 1966, Atlas of Peculiar Galaxies, Pasadena 
Barnes, J. & Hut, P. 1986, Nature, 324, 446 
Bender, R., Burstein, D. & Faber, S.M. 1992, ApJ, 399, 462 
Borne, K.D. & Hoessel, J.D. 1988, ApJ, 330, 51 
Burstein, D. & Heiles, C. 1982, AJ, 87, 1165 
Capelato, H.V., de Carvalho, R.R. & Carlberg, R.G. 1995, 
ApJ, 451, 525 

Ciotti, L., Lanzoni, B. & Renzini, A. 1996, MNRAS, 282, 1 

Davoust, E. & Considere, S. 1995, A&AS, 110, 19 

de Carvalho, R.R. & Djorgovski, S. 1992, ApJ, 389, L49 

de la Ros a I.G., de Carvalh o, R.R. & Zepf, S.E. 2001, AJ, in 

press ( [a,stro-ph/0104324| ) 
Djorgovski, S. & Davis, M. 1987, ApJ, 313, 59 
Dressier, A., Linden-Bell, D., Burstein, D., et al. 1987, ApJ, 

313, 42 

Hernquist, L. 1990, ApJ, 356, 359 

Hernquist, L. 1993, ApJ, 86, 389 

Hickson, P. 1982, ApJ, 255, 382 

Hjorth, J. & Madsen, J. 1995, ApJ, 445, 55 

Karachentsev, I.D. 1972, Catalogue of isolated pairs of galaxies 

in the northern hemisphere, Soobshch. Spets. Astrof. Obs 

7, 1 



8 



E.Evstigneeva, V.Reshetnikov, N.Sotnikova: Fundamental Plane 



Keel, W.C. & Wu, W. 1995, AJ, 110, 129 

Levine, S.E. & Aguilar, L.A. 1996, MNRAS, 280, 13 

McElroy, D.B. 1995, ApJS, 100, 105 

Pahre, M.A., de Carvalho, R.R. & Djorgovski, S.G. 1998, AJ, 
116, 1606 

Shier, L.M. & Fischer, J. 1998, ApJ, 497, 163 
Simien, F. & Prugniel, Ph. 1997a, A&AS, 122, 521 
Simien F. & Prugniel, Ph. 1997b, A&AS 1126, 15 
Teuben, P. 1995, ASP Conf. Ser., 77, 398 
Turner, E.L. 1976, ApJ, 208, 20 

van Dokkum, P.G., Franx, M., Fabricant, D. et al. 1999, ApJ, 
520, L95 

Vorontsov-Velyaminov, B.A. 1959, Atlas and catalogue of in- 
teracting galaxies, Sternberg Institute, Moscow State Uni- 
versity 

Zepf, S.E. & Whitmore, B.C. 1993, ApJ, 418, 72 



